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ABSTRACT: NLR family pyrin domain containing 3 (NLRP3) is a key component of the inflammasome, whose assembly is a crucial 
part of the innate immune response. The aim of the present study was to evaluate the association between exon 3 polymorphisms of 
NLRP3 and the susceptibility to digestive disorders in rabbits. In total, five coding single-nucleotide polymorphisms (cSNPs) were 
identified; all of which are synonymous. Among them, c.456 C> and c.594 G> were further genotyped for association analysis based on 
case-control design (n =162 vs n =102). Meanwhile, growing rabbits were experimentally induced to digestive disorders by feeding a 
fiber-deficient diet, subsequently they were subjected to mRNA expression analysis. Association analysis revealed that haplotype HI 
(the two cSNPs: GT) played a potential protective role against digestive disorders (p<0.001). The expression of NLRP3 in the group 
HlHXi (HlHXi is composed of HlHl, H1H3 and H1H4) was the lowest among four groups which were classified by different t5'pes of 
diplotypes. Those results suggested that the NLRP3 gene was significantiy associated with susceptibility to digestive disorders in rabbit. 
(Key Words: Digestive Disorders, NLRP3, SNP, Rabbit) 



INTRODUCTION 

Pyrin domain containing 3 (NLRP3), coding NLRP3 
protein, is a member of the NALP subfamily in the 
CATERPILLER family (Ting et al., 2006). Several 
members of the cytosolic NLR family have been identified 
as the key regulators of cytokine production (Kanneganti et 
al., 2007). NLRP3 interacts with the adaptor protein 
Apoptosis-associated speck-like protein to activate 
caspase-1 in inflammasomes, which are protein complexes 
responsible for the maturation and secretion of interleukin- 
ip (IL-ip) and IL-1 (Agostini et al., 2004; Dowds et al., 
2004; Lamkanfi et al., 2007). NACHT domain (domain 
present in NAIP, CIITA, HET-E and TPl), also known as 
NBS (nucleotide-binding site domain) coded by exon 3 of 
the rabbit NLRP3 gene, plays a pivotal role in molecular 
recognition of host innate immunity (Royet and Reichhart, 
2003). The variations and functions of NLRP3 gene, 
especially the exon 3 of the NLRP3 gene of mouse and 
human, have been widely investigated (Bauer et al., 2010; 
Roberts et al., 2010). However, the association between 
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genetic polymorphisms and expression of NLRP3 and in 
domestic rabbit digestive disorders has not been reported. 

In recent decades, the incidence of disease has increased 
along with the dramatic genetic improvement of rabbit 
fertility. Among them, digestive disorders are one of the 
major common diseases (Rosell, 2003). Resistance to 
digestive disorders has been proposed to be genetically 
determined in experimental populations (Eady et al., 2007; 
Garreau et al., 2008). In this study, exon 3 of the rabbit 
NLRP3 gene was resequenced for mutational analysis; 
subsequently the association between the NLRP3 mutation 
and the susceptibility to digestive disorders was estimated 
in a control-case population. Meantime, we also 
experimentally induced digestive disorders by feeding a 
fiber-deficient diet to growing rabbits, which were then 
studied for an association between mRNA expression and 
NLRP3 gene genotype. 

MATERIAL AND METHODS 

Ethics statement 

Animal care and tissue collection procedures involved 
in the present study were approved by the Institutional 
Animal Care and Use Committee in College of Animal 
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Science and Technology, Sichuan Agricultural University, 
Sichuan, China (DKY-B20090908). 

Recruitment of control-case group 

In the experimental rabbit farm of Sichuan Agricultural 
University, the recruitment process of case and control 
rabbits has been specified in our previous study (Zhang et 
al., 2011). In brief, after being weaned at 28 d of age. New 
Zealand White rabbits were fed with pelleted food (16% 
protein, 10.8 MJ/kg) until 84 d of age. The food was 
restricted to approximately 80% of average ad libitum 
intake and water ad libitum. During the experimental period 
from the 28th to 84th d, the dead individuals were clinically 
autopsied, and only symptoms of watery content, dilatation, 
impaction, congestion and mucus in the intestinal tract were 
finally recruited into the case group. In the same rabbit 
population, the healthy rabbits without clinical records of 
any disease were collected into the control group. A total of 
264 New Zealand White (NZW) rabbit ear tissues were 
collected from 162 cases and 102 controls in order to 
determine that a genetic relationship within three 
generations didn't exist in the control-case group. 

Experimentally induced digestive disorders 

Fiber deficient diet as one of the most important risk 
factors can induce the incidence of digestive disorders in 
growing rabbit (Gidenne et al., 2000; Bennegadi et al., 
2001). Thus in this study, rabbits were artificially induced 
to a status of digestive disorders by feeding the fiber- 
deficiency diet (CF = 9%, in contrast to the standard diet of 
CF = 15%). Briefly, a total of 60 healthy growing rabbits at 
49 d of age were randomly selected. Following a 7-d dietary 
transition period, all rabbits were fed with the fiber- 
deficient diet (CF = 9%). 

Animals were carefully observed twice a day with 
precise records for all clinical signs of digestive disorders, 
such as reduced feed intake, diarrhea, constipation (caecal 
impaction), and presence of mucus in excreta. On the 70th 
day, all rabbits were slaughtered by intravenous 
administration of sodium pentobarbital and subsequently 



subjected to histopathological examinations and an autopsy. 
In the end, 41 rabbits were successfully sampled and 
classified into three groups according to different severity 
status of digestive disorders based on chnical signs and 
gastrointestinal tract symptoms. The individuals that could 
not be clearly classified into any group were excluded. In 
brief, there were 8 in the healthy group, 9 in the mild 
inflammation group and 24 the in severe inflammation 
group. Colons were collected for RNA isolation and 
immediately frozen in liquid nitrogen and stored at -80°C. 
Ear tissues of 41 rabbits were also collected for DNA 
isolation. 

DNA and RNA isolation 

The ear tissues from induced digestive disorders group 
and control-case group were used for genomic DNA 
extraction by AxyPrep Genomic DNA Miniprep Kit 
(Axygen, USA). The frozen colon tissues from 
experimentally induced digestive disorders group were used 
for RNA extraction by RNAiso Pure RNA Isolation Kit 
(TaKaRa, China). A NanoVue Plus (GE, American) was 
used to assess the concentration and quality of DNA and 
RNA on the basis of absorbance of UV light at 260 (A260) 
and 280 nm (A280). The value range from 1.6 to 1.8 means 
the qualification for next study. 

Screening mutation site in NLRP3 

The genetic polymorphisms in exon 3 of NLRP3 gene 
were scanned in 24 rabbits from control-case group by 
Sanger direct sequencing method. According to the 
reference sequence of rabbit NLRP3 gene (GenBank ID: 
NW_003159516), 2 primer pairs were used for mutation 
screening (Table 1). The PCR reaction condition was 
performed as follows: one denaturation cycle at 94 °C for 5 
min, followed by 39 cycles at 94°C for 30 s, 62 to 62.2°C 
for 30 s and 72°C for 65 s, then an extension cycle at 72°C 
for 10 min. The 30 (il reaction volume included 15 |li1 
2xTaq PCR MasterMix (Aidlab, China), 1.2 |il of each 
primer (10 pmol/^1), 3 nl DNA template (20 ng/|al), 9.6 |al 



Table 1. Inlbrnuuion of Ihc primers used lor PCR, HRM analys: 


is and R'f-PCR 






Primers 


Primer F (5'-3') 


Primer R (5'-3') 


Amplicon size 
(bp) 


Annealing 
temperature 
(°C) 


SNP scanning 










Seq 1 


ACAGGAAGCAGGTTGAAGAAAGG 


GGTGGTGATGAGCAGTGAGGC 


871 


62 


Seq2 


CCACCACCTCCAAGACCAC 


CAGCAGACGCTTACTCTACCC 


922 


62.2 


Genotyping 










HRM456 


GCTGTGCCCGCAGTCCGATT 


AGGCGAAGCCGTGTGTAGCG 


151 


60.5 


HRM594 


GCGGGAGCAGGAGTTGGTGG 


GCTGCCCTGAAGCCCAGTCC 


197 


63.2 


Quantification 










HPRT 


CGAGGACTTGGAAAGGGTG 


CAGCAGGTCAGCAAAGAACT 


151 


63 


NLRP3 


TGAGGAAGAGGACACGGGACG 


GAGCCTGGTGGACCTGATTGC 


92 


58.8 
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ddH20. PCR products were sequenced as previously 
described (Zhang et al., 2011). 

Genotyping using high-resolution melting (HRM) 

Primer pairs of HRM456 F/R and HRM594 F/R (Table 

1) were designed to amplify the small fragment containing 
the target SNPs of c.456 C>G and c.594 G>T, which was 
further subjected to HRM analysis in control-case group. 
PCR reactions were performed on Bio-Rad CFX96 real- 
time PCR detection system (Bio-Rad, Inc., Hercules, CA, 
USA). All samples were amplified in duplicate, and each 
run contained a non-template control (NTC) and three types 
of individuals as the standard including sequencing-verified 
reference sequence-type, heterozygous and homozygous 
genotypes. The PCR reaction used SsoFast EvaGreen 
Supermix (Bio-Rad, CA, USA), the volume and condition 
were similar as previously described (Zhang et al., 2011). 
HRM curve data was analyzed using the manufacturer's 
software. 

Gene mRNA expression in colon tissue using real-time 
RT-PCR 

RNA was used to synthesize cDNA applying the 
PrimeScript RT reagent kit (TaKaRa, China). Two pair 
primers (Table 1) were designed to amplify the short 
fragment including reference gene hypoxanthine 
phosphoribosyltransferase (HPRT, GenBank ID: 
EF062857.1) and the target gene NLRP3 (GenBank ID: 
NW_003159516). PCR reactions were performed on Bio- 
Rad CFX96 real-time PCR detection system (Bio-Rad, Inc., 
Hercules, CA, USA). The total reaction volume of 10 
|Lil/well included 5 |li1 of iQ SYBR Green Super mix (Bio- 
Rad, CA, USA), 0.4 |il of each primer (10 pmol/|il), I |il 
cDNA template (20 ng/^il), and 3.2 ddHjO. Reaction 
condition was 98.0°C for 3 min, 39 cycles of 98.0°C for 5 s 
and 56 to 63.0°C for 25 s, subsequently 95.0°C for 10 s, 
then Melt Curve from 65°C to 95°C per increments of 
0.2°C for 10 s to read plate. The standard curves were 
diluted 10-fold gradient from 10"' to 10"' to ensure an 
amplification efficiency of 100%±5%. Then 41 
experimentally induced digestive disorders samples were 
run in triplicate for the three genes, each run containing a 
NTC. 



Data analysis 

The DNA sequences were assembled and analyzed 
using the DNAstar program (DNAS Inc, Madison, WI, 
USA). The relative synonymous codon usage (RSCU) value 
(SNc/Na) was associated with a change in local translation 
elongation rates (Sharp and Li, 1986), as previously 
described to assess the synonymous mutation (Sauna and 
Kimchi-Sarfaty, 2011). The alleles and genotypes frequency 
differences between control-case groups were tested with 
the Fisher's exact test. Odds ratios (ORs) for the allelic 
association were estimated by unconditional logistic 
regression analysis. 

The individual gene expressions were quantified via 
real-time RT-PCR and normalised by the vandasompele 
method (Vandesompele et al., 2002). Gene expression data 
was assessed for normal Gaussian distribution and then 
analyzed by ANOVA test and LSD test. Haplotypes were 
constructed with PHASE version 2.1 with the default 
parameters (Stephens et al., 2001; Zheng et al., 2004). The 
frequency differences of alleles and genotypes between the 
case group and control group were tested with the Fisher's 
exact test. Hardy-Weinberg equilibrium (HWE) and pair- 
wise LD tests were performed in SHEsis web-based 
platform (Li et al., 2009). All statistical analyses were 
performed using SAS v8.0 statistical package (SAS 
Institute, NC, USA). p<0.05 was considered statistically 
significant. 

RESULTS 

Mutation analysis of rabbit NLRP3 gene 

Exon 3 of rabbit NLRP3 gene was detected in 24 
individuals (12 individuals were randomly selected from 
case group and control group, respectively). The sequencing 
result showed that five coding single nucleotide 
polymorphisms (cSNPs) (c.456 C>G c.594 G>T, c.1224 
G>T, C.1248 G>A, c.1627 G>A) were all synonymous 
mutations. In order to represent the effect of synonymous 
mutations on codon usage bias, the RSCU value was 
analyzed (Table 2). ARSCU values in c.456 C>G (-24.7), 
c.594 G>T (-8.2) and c.1627 (-1.4) were negative (ARSCU 
<0), which means the variation introduces a rarer codon 
associated with a slower local rate of translation elongation 



Table 2. Codon usage of synonymous mutations in the coding NLRP3 gene 



Codon change RSCU change . . 

mRNA position ARSCU Amino acid 

From To From To 

456 GCC GCG 342 9^5 ^247 A 

594 GGG GGT 17 8.8 -8.2 G 

1224 CCG CCT 8.7 12.6 3.9 P 

1248 ACG ACA 9.1 11.7 2.6 T 

1627 AGG AGA 10.6 9.2 -1.4 R 



RSCU = Relative synonymous codon usage. 
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compared with the reference type. While the ARSCU values 
in C.1224 G>T (3.9) and c.1248 G>A (2.6) were positive 
(ARSCU>0). Moreover, the higher ARSCU value is more 
suitable for "optimal codon" (Comeron, 2004). Among the 
five mutations, the ARSCU value of c.456 C>G (-24.7) and 
C.594 G>T (-8.2) were obviously higher than others. Thus, 
synonymous mutations c.456 C>G and c.594 G>T were 
chosen for HRM analysis. 

Case-control association of c.456 C>G and c.594 G>T 
mutations 

Total of 264 NZW rabbits including 162 cases and 102 
controls were genotyped using HRM technology in c.456 
C>G and c.594 G>T site. The frequencies in case vs control 
group were significantly different in allele C (c.456 C>G 
48% vs 30%, p<0.01) and allele G (c.594 G>T, 62% vs 
53%, p<0.01). Similarly, the frequencies of GG genotype 
(16.7% vs 50%, p<0.01) and TT genotype (8% vs 24.5%, 
p<0.01) were significantly different. The case-control 
association analysis revealed that both allele G and T 
carried a potential protective role for digestive disorders 
with odds ratio (OR) value of 0.4537 (case vs control, 95% 
confidence interval 0.3132 to 0.6573, p<0.01) and 0.6975 
(case vs control, 95% confidence interval 0.4892 to 0.9944, 
p<0.05) respectively. Under these models, association 



analysis also suggested GG genotype decreased risk of 
digestive disorders with an odds ratio (OR) value of 0.20 
(case vs control, 95% confidence interval 0.11 to 0.35, 
p<0.01) while G/G-G/T genotype increased risk of digestive 
disorders with an odds ratio (OR) value of 3.72 (case vs 
control, 95% confidence interval 1.80 to 7.68, p<0.01) 
(Table 3). 

The haplotype analysis showed that two cSNPs had 
formed four haplotypes among the case and control groups 
composed of 264 individuals. The four haplotypes were 
conveniently defined as allele HI (the two cSNPs: GT), 
allele H2 (the two cSNPs: CG), allele H3 (the two cSNPs: 
GG) and allele H4 (the two cSNPs: CT) to simplify the term 
usage. The case-control association analysis revealed that 
haplotype HI played a potential protective role against 
digestive disorders (Table 4). 

Among the 264 individuals of case and control groups, 
the average frequencies of haplotypes HI, H2, H3 and H4 
were 0.3568, 0.3530, 0.2303 and 0.0599 respectively (Table 
4), which further determined the diplotypes frequencies of 
4.9% (HlHl), 33.7% (H1H2), 13.6% (H1H3), 7.2% 
(H1H4), 2.7% (H2H2), 17.4% (H2H3), 7.6% (H2H4), 
11.0% (H3H3) and 1.9% (H4H4). The HWE test did not 
reject the null hypothesis (p>0.05). Significant deviations 
from Hardy-Weinberg equilibrium (HWE) were observed 



Table 3. NLRP3 alleles and genot5T3es associated with digestive disorders in rabbit 


Models 


Genotypes/aUele 


Cases (%) n = 162 


Controls (%) n = 102 


OR (95% CI) 


p values 




C.456C 


157 (48%) 


61 (30%) 


1.00 


<0.0001 




C.456G 


167 (52%) 


143 (70%) 


0.4537(0.3132-0.6573) 




Co-dominant 


G/G 


27 (16.7%) 


51 (50%) 


1 


<0.0001 




G/C 


113 (69.8%) 


41 (40.2%) 


0.19(0.11-0.35) 






c/c 


22 (13.6%) 


10 (9.8%) 


0.24(0.10-0.58) 




Dominant 


G/G 


27 (16.7%) 


51 (50%) 


1 


<0.0001 




G/C-C/C 


135 (83.3%) 


51 (50%) 


0.20(0.11-0.35) 




Recessive 


G/G-G/C 


140 (86.4%) 


92 (90.2%) 


1 


0.35 




C/C 


22 (13.6%) 


10 (9.8%) 


0.69 (0.31-1.53) 




Over-dominant 


G/G-C/C 


49 (30.2%) 


61 (59.8%) 


1 


<0.0001 




G/C 


113 (69.8%) 


41 (40.2%) 


0.29(0.17-0.49) 




Log-additive 








0.35 (0.23-0.55) 


<0.0001 




C.594G 


200 (0.62) 


108 (0.53) 


1.00 


0.0463 




C.594T 


124 (0.38) 


96 (0.47) 


0.6975(0.4892-0.9944) 




Co-dominant 


G/G 


51 (31.5%) 


51 (31.5%) 


1 


0.0008 




G/T 


98 (60.5%) 


98 (60.5%) 


0.77 (0.44-1.36) 






T/T 


13 (8%) 


13 (8%) 


3.16 (1.41-7.08) 




Dominant 


G/G 


51 (31.5%) 


51 (31.5%) 


1 


0.85 




G/T-T/T 


111 (68.5%) 


111 (68.5%) 


1.05 (0.62-1.80) 




Recessive 


G/G-G/T 


149 (92%) 


149 (92%) 


1 


0.0002 




T/T 


13 (8%) 


13 (8%) 


3.72 (1.80-7.68) 




Over-dominant 


G/G-T/T 


64 (39.5%) 


64 (39.5%) 


1 


0.014 




G/T 


98 (60.5%) 


98 (60.5%) 


0.54 (0.32-0.89) 




Log-additive 








1.51 (1.03-2.23) 


0.033 



OR = Odds ratio, CI = Confidence interval. 
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Table 4. Haplotype association witli digestive disorders in the 
rabbit (n = 264) 



Haplotype 


Cases/total 


OR (95% CI) 


p-value 


HI 


23/81 


1.00 




H2 


57/80 


6.24 (3.15-12.38) 


<0.001 


H3 


144/229 


4.27 (2.46-7.42) 


<0.001 


H4 


100/138 


6.63 (3.60-012.22) 


<0.001 



OR = Odds ratio, CI = Confidence interval. 

only in control group (c.456:p 
0.00036). 



4.53x10'^ c.594:p 



NLRP3 gene mRNA expression in digestive disorders 

Forty-one digestive disorders individuals consisted of 8 
in healthy group (H), 9 in mild digestive disorders group 
(M) and 24 in severe digestive disorders group (S). The 
mRNA expression of NLRP3 gene was detected using real- 
time RT-PCR. The mRNA expression of the NLRP3 gene in 
different states of digestive disorders is exhibited in Figure 
1. The expression level of the NLRP3 gene in group S 
(0.3455+0.0440) was significantly higher than that in group 
M (0.0452±0.0314, p<0.05), and dramatically higher than 
that in group H (0.2256+0.0193 p<0.01). 

NLRP3 gene mRNA expression in different genotypes of 
NLRP3 gene 

The NLRP3 gene was genotyped by HRM analysis in 41 
rabbits which were experimentally induced to digestive 
disorders. The different expression in six genotypes of 
NLRP3 genes are shown in Figure 2. In SNP c.456 C>G the 
NLRP3 gene expression level of CC genotype 
(0.6379+0.1354) was significantly higher than GC genotype 
(0.1526+0.0650) and GG genotype (0.1211+0.1105, p< 
0.01). In SNP C.594 G>T the NLRP3 gene expression level 



0.45 
0.4 

I 0.35 

I 
a. 

S 0.3 

« 

s 

£ 0.25 
a: 

0.2 
0.15 
0.1 

0.05 



p=o.o2es 



s(rg=24) 



Figure 1. The expression of NLRP3 in different severity states of 
digestive disorders. H, M, S and N respectively represent healthy 
group, mild digestive disorders, severe digestive disorders group 
and the number of each group. The histogram represents the 
results as mean±standard error of the mean (SEM). Only the 
significant level (p<0.05) can be marked out. 

of TT genotype (0.0302+0.0191) was significantly higher 
than GG genotype (0.3470+0.0923, p<0.05). 

The mRNA expressions of the NLRP3 gene were 
divided into four groups according to the different 
diplotypes of the 41 rabbits. Group HlHXi was composed 
of HlHl, H1H3 and H1H4. Group H2HX2 was composed 
of H2H2, H2H3 and H2H4. Group HX3HX3 was composed 
of H3H3 and H4H4 (Figure 3). The result showed that the 
NLRP3 expression in group HIHX was the lowest among 
four groups and significantly lower than that in group 
H2HX (p<0.0001). 



0.9 

0.8 

o 0-7 H 

2 0.6 
a. 

2 0.5 

g 0.4 
5 0.3 

^ 0.2-1 

0.1 
0 



P =0.0053 



P=0.0025 



P=0.0423 



CC(N=6) I GC(N=26) | GG(N=9) 
c.456 



TT(N=9) 



GT(N=17) 
c.594 



GG(N=15) 



Figure 2. The expression differences of NLRP3 in different allies for NLRP3 genotypes. N represents the number of each allies. Values 
are presented as mean+standard error of the mean (SEM). Only the significant level (p<0.05) can be marked out. 
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P =0.0032 





H1HX1{N=12) H1H2(N=12) H2HX2(N=14) HX3HX3(N=3) 

Figure 3. The expressions of NLRP3 gene in different diplotypes of rabbits. Group HlHXl is composed of HlHl, H1H3 and H1H4. 
Group H2HX2 is composed of H2H2, H2H3 and H2H4. Group HX3HX3 is composed of H3H3 and H4H4. N represents the number of 
each group. The histogram represents the results as mean±standard error of the mean (SEM). Only the significant level (p<0.05) is 
marked out. 



DISCUSSION 

Digestive disorder is a complex multifactorial disease, 
the pathogenesis of which is still not well understood. 
NLRP3 is a member of the CATERPILLER family of genes 
encoding proteins that comprise a nucleotide -binding 
domain and a leucine-rich repeat domain. Cryopyrin, the 
protein encoded by NLRP3, controls the inflammasome, a 
crucial molecular platform that regulates activation of 
caspase-1 and processing of interleukin (IL)-ip, two key 
mediators of inflammation (Martinon and Tschopp, 2006; 
Mariathasan and Monack, 2007; Petrilli et al., 2007). 
Mutations in NLRP3 are known to be responsible for three 
rare autoinflammatory disorders. The combination of the 
polymorphisms CARDS (ClOX) and NLRP3 (Q705K) was 
shown to be associated with rheumatoid arthritis (Kastbom 
et al., 2008). Men who have both the ClOX and Q705K 
alleles in CARDS and NLRP3, and who express wild-type 
alleles of Nod2 are at an increased risk of developing 
Crohn's disease (CD) (Schoultz et al., 2009). CD is an 
idiopathic inflammatory bowel disease in humans 
comparable in many features to bovine paratuberculosis 
involving an abnormal mucosal immune response. 
Cummings et al. (2010) found four SNPs in NLRP3 are in 
high linkage disequilibrium (LD) with each other. 
Rs4925648 and rsl0925019 are the most strongly 
associated with CD susceptibility. Therefore, NLRP3 is 
considered as one of the most plausible candidate genes 
associated with digestive disorders in rabbits. 

The degeneracy of the genetic code means that many 
mutations in coding sequences, especially at the third base 
of codons, do not affect protein sequence and are therefore 
considered silent. However, there is increasing evidence 



that they could still have effects on transcription, splicing, 
mRNA transport or translation, any of which could alter the 
phenotype (Goymer, 2007) and cause diseases. In the 
present study, we, for the first time, found the NLRP3 
polymorphisms was significantly correlated with the 
genetic resistance of rabbits digestive disorders, which 
expands our knowledge about the important role of NLRP3 
polymorphisms in determining the genetic resistance to 
infectious diseases in domestic animal based on case- 
control study. Although these coding SNPs of c.456 C>G 
and C.594 G>T did not result in amino acid change, the 
association analysis also revealed that both allele c.456 G 
and allele c.594 T were significantly associated with 
decreased risk of developing digestive disorders in rabbits 
(p<0.05). In addition, the haplotype analysis also supported 
this point of view. However, the potential biological effects 
or causative variants linked to this synonymous SNP should 
be further investigated. 

Among the 162 case subjects recruited from the general 
commercial population, we guaranteed their homogeneity 
by performing the same selection criteria according to the 
clinical manifestations and pathological symptoms, which 
has been described in detail in our former study (Zhang et 
al., 2011). However, we could not absolutely exclude the 
possibility that these case individuals could be further 
divided into different subtypes of digestive disorders by 
subjecting them to intestinal microbiota analysis with 
cultivation-based and nucleotide sequence-based methods 
(Andoh et al., 2009). Furthermore, we were also unable to 
confirm whether an individual death was caused by primary 
digestive disorders due to the underdeveloped clinical 
diagnosis in farm animals. This potential heterogeneity 
among case subjects would more or less reduce the 
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detection power for our case-control association study. Even 
so, mRNA expression analysis also revealed the practical 
roles of NLRP3 involved in the development of digestive 
disorders, which was experimentally induced by feeding a 
fiber-deficient diet. In rabbits, the dietary fiber level has 
been considered as an important environmental factor 
associated with digestive disorders. However, the direct 
influence of dietary fibre component on NLRP3 mRNA 
expression has not been reported. 

Strong deviations from HWE are believed to most likely 
associate with genotyping errors (Ott, 2004). However, the 
significant deviations from HWE for c.456 C>G and c.594 
G>T variants observed in the present study could not be 
explained by potential genotyping errors because we 
performed strict quality control for HRM analysis (see 
Materials and Methods for details). Because there is 
selective pressure (Digestive disorders) in the case group, it 
is reasonable that non-random mating would be the main 
factor accounting for the failure in fitting assumption of 
HWE (Liu and Hu, 2010). 

In our study, the expression of NLRP3 significantly up- 
regulated following the aggravation of intestinal 
inflammation (Figure 1), this result is similar to an earlier 
study that found the expression of the NLRP3 gene is low in 
colons of healthy mice but elevated in colons of mice with 
colitis (Villani et al., 2008). Dupaul-Chicoine et al. (2010) 
also found that loss of NLRP3 resulted in more severe DSS 
colitis and Hirota et al. (2011) found that Nlrp3"'"mice were 
more susceptible to experimental colitis. At the same time, 
we found that the NLRP3 expression level of genotype GG 
(c.456) and TT (c.594) in severe digestive disorders status 
were significantly and respectively lower than that in 
genotype CC (c.456) and GG (c.594) (Figure 2). The 
expression of NLRP3 in the group HlHXi was the lowest in 
the four groups and highly significantly lower than that of 
the group H2HX2 (p<0.0001) (Figure 3). Of course, the 
potential biological effects on cytokines expression or 
causative variants linked to this synonymous SNP should be 
further investigated. 

CONCLUSION 

Both the genetic association study and the gene 
expression analysis suggest that the NLRP3 gene plays an 
important role in the susceptibility of digestive disorders. 
Thus, NLRP3 has a potential as a candidate gene for 
digestive disorders in rabbit. However, further research is 
required to confirm this conclusion and explore the 
potential biological implications. 
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